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Transition of the ground state of a classical $ 4 theory in 2 + 1 dimensions is 
studied from a metastable state into the stable equilibrium. The transition occurs 
in the broken Z2 symmetry phase and is triggered by a vanishingly small amplitude 
homogeneous external field h. A phenomenological theory is proposed in form of 
an effective equation of the order parameter which quantitatively accounts for the 
decay of the false vacuum. The large amplitude transition of the order parameter 
between the two minima displays characteristics reflecting dynamical aspects of 
the Maxwell construction. 



The range of interest of the irreversible decay of a metastable vacuum 
state of finite energy density covers effects from cosmological phase transitions 
to instabilities observed im±he mixed phase of first order phase transitions of 
condensed matter systems El. 

Whether the relevant mechanism for a first order phase transition is the 
formation of bubbles of the new phase, as described by thermal nucleation 
theory, or the gradual change of a large region of the sample, due to small 
amplitude spinodal instabilities described by spinodal decomposition is also 
an intriguing question in heavy ion physics where^the actual expansion rate of 
the plasma may favour one or the other scenarioB. 

The conventional treatment of the decay of metastable states is based on 
the nucleation theory but concurrent small amplitude spinodal instabilities 
are also present in the system. They are responsible for the flattening of the 



static effective potential (Maxwell cut) U. The clarification of their role in the 
metastable — > stable transition is the main theme of the present investigation. 

1 The model and the time-history of the order parameter (OP) 

We study numerically the dynamics of a classical $ theory in 2 + 1 dimensions 
governed by the discretized field equation of motion + a t ) + $n(t — a t ) — 

2* A (t) + a?(-<S>n +*%-h)-$ Ei(*A+j(*) + ~ 2 *™ (*)) = 0, with 

initial conditions : <&n(t = 0) = 0, $&(t = 0) = $0 + £n$i where £ft is 

^Presented at Strong and Electroweak Matter (SEWM2000), Marseille, France, 14-17 June, 
2000. Based on an article appeared in Phys. Rev. D62 (2000) 085013 written in collabora- 
tion with Szabolcs Borsanyi, Andras Patkos, Janos Polonyi. 
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Figure 1: Main qualitative features of the <3>(t)-trajectory for a 512 X 512 lattice and the 
external field h = -0.04/v^. 

an evenly distributed white noise in the range (—1/2, 1/2) (all quantities are 
expressed in units of mass). The corresponding initial kinetic power spectrum 
is E k (k) ~ w 2 (k) = -1 + 4(sin 2 (fc K a/2) + sm 2 (k y a/2))/a 2 . We have chosen 
$0 = 0.815 and <&i = 4/ y/6 which corresponds to a temperature value Tj = 0.57 
assuring that the system is in the broken symmetry phase. Our goal is to 
describe the evolution of the system only in terms of an effective equation for 
the OP $(*) = £ 

Fig. [j] shows a typical OP-history together with the later time history of 
the OP mean square (MS)-fluctuation (( < & 2 ) — ($) 2 ) and of its third moment 
(((<& — (<I>)) 3 )). In general one can distinguish five qualitatively distinct parts 
of the OP-trajectory that starts with large amplitude damped oscillations cor- 
responding to the excitation of resonating modes, followed by a rather slow re- 
laxation to a metastable state characterised by (<&) ~ 0.72. The (quasi) thermal 
motion in the metastable state is followed by the metastable — > stable transi- 
tion induced by the external field h, during which we can observe characteristic 
variations of the second and third moments. Quantitative interpretation of this 
variation will be presented in the following section. The last portion of the tra- 
jectory represents thermal motion in the true ground state. 

2 The effective OP-theory 

For the description of (quasi)thcrmal motion near a (meta)stable point we 
assume the validity of the ergodicity hypothesis for a system which consist of 
a single degree of freedom, the OP of the system. We describe its local time 
evolution by an effective Newton type equation: i&(t) +r) (t) + f (&) = £(t). 
ry(<i>), /(<&), C(i) are obtained by a fitting procedure attesting in this way also 
the presence of a term violating time-reversal invariance ?7($). C(0 is the 
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Figure 2: Comparison of the force measured on a 40 X 40 lattice with te one calculated from 
a one-loop and a two-loop effective potential. On the right hand one can see the bending 
down of the force prior the transition time for a 100 X 100 lattice and h = —0.04/^/6. 



"error" of the best global fit to the homogeneous equation at time t. The 
force felt by the OP, /(<&) agrees with the force derived from the equilibrium 
two-loop effective potential as shown in Fig. |[ To probe this agreement in 
a relatively wide region one has to measure the force by shifting the center 
of motion to different values of $ by applying appropriate h fields to the 
system. The coefficient rj is well-defined and positive, but its value depends 
quite substantially on the time resolution. The analysis of the relationship of 
time averaging to the nonzero value of r] is left for future investigations. 

The thermalization "history" of the system is shown in Fig. g[ The dis- 
played temperature variations correspond to the kinetic energies of the soft and 
hard modes. During the time evolution they approach each other, manifesting 
relaxation to a local thermal equilibrium. Because of our "white noise" ini- 
tial condition the equilibration process is characterized by an energy transfer 
towards the low k modes. 

The fluctuation moments depicted in Fig. [l] along with the measured force 
shown on the right of Fig. || tell about how the transition proceeds. The in- 
creased values of the moments indicate the enhanced importance of the soft 
interactions. Preceding directly the transition towards the direction of neg- 
ative <!> values, the fitted force bends down and its average becomes a small 
positive constant. The vanishing of the force implies the flatness of the effective 
potential along the motion of the OP (the mode k = in momentum space) 
indicating the dynamical realization of the Maxwell-cut. 

Taking into account the existence of a mixed phase during the transition 
period we can construct a model that reproduces exactly the shape of the mo- 
ments and shows the vanishing of the force felt by OP. Upon space averaging 
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Figure 3: The time evolution of the kinetic energy content of the |k| > 2.5 and |k| < 2.5 
regions averaged over the corresponding |k|-intervals. 



the microscopic field equation of motion of section 1 we find the following equa- 
tion of motion for the OP: = $ - $ + $ 3 + 3<^' $ -h = ® + dV mst /d$ 
where the symbol tp nV means the space average of ip n , (Tp v = 0,$(i,x) = 
+ (p(x,t)). The instant potential Vi ns t, which is a fluctuating quantity, 
contains as a deterministic piece the sum of the tree-level potential and of the 
slowly varying part of the second and third moments <p nV (t) — ((p nV )det+Cn(t), 
n = 2, 3. We assume that the space can be splitted into sharp domains (ne- 
glecting the thickness of the walls in between), where the field is the sum 
of the constant background values <&o± and the fluctuations tp± around it, 
*J>±(x, t) = <£>o± + <p± (x, t). Based on the smooth evolution of the temperature 
as displayed in Fig. [| we assume also local equilibrium in both phases. 

The actual value of the OP is determined by the surface ratio p(t) occupied 
by the stable phase: <&(i) = p(t)$>o- + (1 — p(t))^o+- Simple calculation yields 

-3<S>(t)(¥ V )det(t)-<5> 3 (t). 
If one takes the values of <I>o± ; ¥>± from the respective equilibria a quite 
accurate description of the shape of the two fluctuation moments arises in the 
whole transition region and its close neighbourhood using the measured to 
parametrize their i-dependence. Substituting the expressions of the moments 
into the equation of the OP one finds for the deterministic part of the force, 

/ (*) = $0+ - $0 _ $ - * + g, 0+ - $0 _ h. The average of 

the equations of motion in the respective equilibria, ($±(x, t)) — $o± — h = 
implies the vanishing of the deterministic force, when exploiting the equality of 
averaging over the volume and the statistical ensemble. The equation of motion 
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for OP goes into $(f) + ( 3 (t) + 3( 2 (t)$(t) = which reflects the dynamical 
realization of the Maxwell construction when assuming the validity of local 
equilibrium in both phases. 



3 The nucleation picture 

The statistics of the release time, the time necessary for the system to escape 
from the metastable region by nucleating a growing bubble of the stable state, 
is of the form P(t) ~ exp(-ir(/i)L 2 ) where T(h) = Aexp(-S 2 (h)/T) is the 
nucleation rate and L is the lattice sizeD. We extract the exponent S 2 by fitting 
the measured rate to the expression of T(h) assuming the /i-independence of 
A. Comparing with the bounce action of the nucleation theory we obtain: 
Sthmwaii = ^§-hT ~ 155 meaeured . Using the T ^ effective potential instead of 
the bare one one can achieve a much better agreement, but still a factor of 2 
discrepancy remains with the mean field approach. 



4 Nucleation vs. spinodal instabilities 

Our results offer a "dualistic" resolution of the competition between the nucle- 
ation and the spinodal phase separation mechanisms in establishing the true 
equilibrium. We find that the statistical features of the decay of the false 
vacuum agree with the results of thermal nucleation. Alternatively, the effec- 
tive OP-theory displays the presence of soft modes and produces dynamically a 
Maxwell cut when the time dependence of the transition trajectory is described 
in the effective OP theory. During the transition the OP travels trough a nar- 
row but flat valley around the k = mode along which we expect the effective 
potential to be flat. This flatness is reflected also in the decrease of the kink- 
like action S measure d relative to Sthinwaii ■ The larger is the system the smaller 
is the external field which is able to produce the instability. 
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